Scaling properties of delay times in one-dimensional random media 
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The scaling properties of the inverse moments of Wigner delay times are investigated in finite 
one-dimensional (ID) random media with one channel attached to the boundary of the sample. We 
find that they follow a simple scaling law which is independent of the microscopic details of the 
random potential. Our theoretical considerations are confirmed numerically for systems as diverse 
as ID disordered wires and optical lattices to microwave waveguides with correlated scatterers. 

PACS numbers: 72.15.Rn, 72.20.Dp, 73.23.-b 



The study of the statistical properties of Wigner de- 
lay times has been a subject of intense research ac- 
tivity 0, 1 1 I, S 0, 0, i, ©, EI EL El El Ei. 

The Wigner delay time is defined as the energy deriva- 
tive of the total phase of the scattering matrix S i.e., 
tw = — ihd In dot S/dE, and can be interpreted as a 
time delay in propagation of the peak of the wave packet 
due to scattering interference, in comparison to a free 
wave packet propagation. Although most of the contem- 
porary activity has been focused in understanding the 
statistical properties of delay times within chaotic meso- 
scopic systems 0, E[, recently the interest has shifted 
towards random scattering media exhibiting Anderson 
localization @, [J 0, 0, [J 0, 0, EI E3 including the 
most difficult case of the Anderson Metal-Insulator Tran- 
sition (MIT) JlE EI El, [IE El- On the experimen- 



tal side d, [i, IfJ, the statistics of scattering phases 
and delay times have been measured in microwave ex- 
periments with quasi-one-dimensional random samples, 
while on the theoretical side the main effort has been 
to connect the statistical properties of delay times with 



that of eigenfunctions [ill . El, EH- Establishing 



such a relation may open new exciting opportunities 
for measuring the statistic al prop erties of eigenfunctions 
0, EE E3, EE El, H ED, S El, Hi via the experimen- 
tally accessible delay times. 

Specifically, usin g th e powerful Non-Linear a Model 
(NLerM) technique [111, [l2j an exact relation was found 
linking the probability distribution of eigenfunction com- 
ponents within a random medium to the distribution of 
Wigner delay times in the same sample of length L, with 
one channel attached at its bulk. This relation is exact 
on the level of the NLaM and valid independent of the 
system size L (i.e. irrespective if we take the thermody- 
namic limit L — > oo or keep L finite). 

However, one has to question the validity of mapping 
a particular microscopic model of a disordered system 
onto the NLerM. More specifically, this mapping is ap- 
proximately correct in the case of weak disorder and 
breaks down totally for strong disorder. Another strict 
requirement is that the underlying geometry allows for 
a diffusive process - this certainly is not the case for 



strictly one-dimensional (ID) random media. Finally, 
NLerM calculations pre-assume that the disorder poten- 
tial is white-noise, thus excluding the emerging family of 
disordered systems with imp rinted correlations in their 
potential 0, EI E3, El, EI H EH, E2 . The above re- 



strictions cast reasonable doubts on the validity of NLerM 
predictions, as far as realistic systems are concerned, and 
call for testing by means of a dedicated experiment or 
computer simulation. 

It is the purpose of this paper to investigate the scaling 
properties of moments of delay times and compare them 
with the ones found for wavefunctions in cases where the 
conditions for NLerM applicability are violated. To this 
end we will study various microscopic systems: (a) a ID 
disordered electronic system (modeled by an Anderson 
Hamiltonian) , (b) a microwave system with long-range 
correlated scatterers inside a waveguide (modeled by a 
Kronig-Penney model), and (c) cold-atoms in a disor- 
dered optical lattice (modeled again by a Kronig-Penney 
model with binary distribution). In all cases, we find that 
the inverse moments of delay times T£ in a disordered 
sample of length L follow a simple scaling law which is 
independent of the microscopic properties belonging to 
the underlying physical system. Specifically we find that 



/?_, = /(A_ g ) ; /3_,= 



_ 'rcf 



" {rZ q Y 



_ 'rcf 



(1) 



where q takes positive values and (tl) represents the av- 
erage (or typical) delay time over disorder realizations. 
The variable represents the delay time of the L — > oo 
sample with the same disordered potential. The variable 
r re f is the delay time of a " reference" sample, correspond- 
ing to an "infinite" localization length set-up, with length 
L. The former quantity incorporates the microscopic in- 
formation of the system (i.e. disorder potential), whereas 
T re f only depends on the information of the finite sample 
length L, as well as the dimensionality and the energy 
E at which the scattering experiment is performed. Our 
numerical analysis indicates that the scaling law, Eq.([T]), 
can take the model-independent form 
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where e is the disorder strength of the random potential. 
In fact, our numerical data, suggest that Eq.([2|) is exact 
only for q = 1 while for higher q— values small deviations 
from the linear behaviour can be detected. 

We point out that a similar relation to Eq.@ was 
found for the scaling properties of wavefunction moments 
within a closed disordered sample [11 G1, E3] ■ Thc 

corre- 
sponding expression involves the q' = q + 1 wavefunction 
moment and reads 
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where l ( p = L(P(«')/P ] ^ ) ) 1 /( 1 -?') are the various in- 
formation lengths of a sample with length L, P^ q ' = 
IV'nl 2 * with eigenfunction components ip n , loo(e,E) 
is the localization length of the infinite sample with the 
same disordered strength, and P' f ~ L with a pre-factor 
defined by the reference geometry (ID periodic lattice in 
the cases studied here). For the special case q' = 1 the 
corresponding information length is equal to the entropic 
length defined by = e/2exp(- £^ =1 \ipn\ 2 In |Vn| 2 )- 

We initiate our analysis by recalling the notion of delay 
times as originally proposed by Wigner [3I [3(| . This is 
the time that a reflected particle is delayed due to inter- 
action with a scattering region. Now we recall that the 

(2) 

q = 2 information length, l L (e, E), (associated with thc 
inverse participation ratio) measures the "penetration" 
(localization) length inside a disordered sample before 
the particle is reflected back (we are considering here the 
one channel scattering set-up). The corresponding delay 
time due to the scattering from the disordered sample is 
then given by tl = 21 L fv, where v is the group velocity 
of the wavepacket centered around energy E. Using this 
argument and substituting it for in Eq.Q, we obtain 
Eq.@ for q = 1. In fact, our numerical data (see below) 
indicate that Eq. @ describes to a good approximation 
higher g-moments as well. 

Below we report our numerical results for various mi- 
croscopic models which support the scaling of Eq.([T|). 
Although our presentation focuses on the first moment 
q = 1, we have found that higher moments follow the 
scaling law, Eq.([T]), equally well. 

ID Disordered Electronic System - The standard 
model that describes a one-dimensional disordered elec- 
tronic sample is the tight-binding equation 



ip n+ i + il) n -i = (E(k) - V n ) ip n ; 



1,2,. ...L (4) 



where k is the incident wavenumber and ip n is the wave- 
function amplitude at the n th site. The on-site potential 
V n for 1 < n < L is independently and identically dis- 
tributed with a box probability distribution, i.e. the V n 
are uniformly distributed on the interval [— e/2, e/2]. 

We open the sample by attaching one channel to the 
first site n = 1. Thc Wigner delay time of a sample of 



length n + 1 is then evaluated with the use of the Hamil- 
tonian map approach Q through the following iteration 
relations 

r n+ i = Gn 1 (r„ + + 1+ [ tan ^ n 2 fc ) +j4n ]!i l°nk' ( 5 ) 
G n = 1 + A n sin [2(0„ - k)] + A 2 n cos 2 (0 n - k), 

where A n = 4^ and the scattering phase is given by 

tan(</>„ + i) = tan(0„ - k) + A n . (6) 

In Fig. [1] we report the delay times for the Anderson 
model, Eq.((4]). The data are averaged over an ensem- 
ble of 10 4 realizations of the random potential and are 
plotted according to the scaling, Eq.Jl}. The value of 
( T oo 9 ) was calculated for a sample of length L = 10 7 
and its convergence was checked by increasing the sys- 
tem size by an additional order i.e. L = 10 8 . The 
scaled data - for various Us and disordered strengths 
e - falls on a single curve, confirming the validity of the 
theoretical prediction, Eq.([T]). Within the same figure 
we also report the corresponding scaled entropic lengths 
(see solid red symbols) / L versus the localization pa- 
rameter A = 2l 00 {e,E)/L, in order to compare with thc 
scaling law that dictates the delay times. The agreement 
between information lengths and delay times is evident, 
thus confirming that these two quantities are directly re- 
lated. In thc inset of Fig. [T] we also report our numerical 
results for the second moment, i.e. q = 2. A nice agree- 
ment with the theoretical expectation, Eq.([T]), is again 
quite evident. 

Microwaves propagating in a ID waveguide - The cre- 
ation of frequency pass/stop bands separated by mobil- 
ity edges and their manipulation by imposing appropri- 
ate correlations in the disordered potential [26| |27|, [28|, 
HSl il| have recently gained considerable research in- 
terest due to their immediate technolo gica l applications. 
One prominent theoretical suggestion [261 ] was based on 
the introduction of long-range correlations in the on-site 
disordered potential. Thc theoretical predictions were 
further supported by subsequent experimental microwave 
measurements [27j |. carried out in a single-mode waveg- 
uide with correlated scatterers realized by screws extend- 
ing from a waveguide wall. By arranging the lengths of 
the screws according to a predefined sequence, correlated 
scattering arrangements could be realized leading to pre- 
defined mobility edges. If the screws are approximated 
by delta scatterers, the propagation of a single mode in 
the waveguide can be described by the wave equation for 
thc Kronig-Penney model 

00 

il>"(z) + Ei/>(z)= e n ip(z n )S(z-nd), (7) 

n— — 00 

where d is the distance between nearby scatterers, ip is 
the electric field of the TE mode, and the energy is given 
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FIG. 1: (Color online) Scaled inverse delay times,Eq.([TJl, for 
the Anderson model. Various symbols correspond to differ- 
ent disordered potentials e £ {0.1,0.5,1,5,10} and \E(k — 
\ZtF)\ < 1. Blue hollow symbols denote delay time data for 
q — 1. For comparison, red solid symbols denote q' = 2 in- 
formation length data, i.e. '/L versus A = 2l 00 {e, E) / L. 
The dashed line is the result of the best fit from Eqs. (|12ll3fl . 
Inset: Same as in the main figure but now q = 2 for delay 
times and q — 3 for information lengths. 



by E = k 2 . We can rewrite the above equation in the 
discrete form for ip n = ip{z n = nd) 

i/Wi + V-Vi-i = [2cos(fcc?) — U n ■ kdsm(kd)] ip n (8) 

One can split the potential U n into a mean e and a fluc- 
tuating term e„, U n = e + e n . Eq. (JSJ) is then equiva- 
lent to the tight-binding equation ([3]), with energy E — > 
2 cos k + ke sin fc and random potential V ra — ► fce n sin fc. 
By choosing the on-site potential as 25|, 2(| 



rn— — oo 



Cm ' C 



n+m 



(9) 



where Cn+m is a random variable, uniformly distributed 
within the interval (0, 1], and 



f(M2-m) 3/2 , 



[sm(2m/x 2 ) - sin(2m/xi) 



TO = 

( 

TO ^ 



(10) 



with /ii = 0.27T, /i2 = 0.47T, and e = —0.1 it was ar- 
gued that mobility edges can be tailored at wavenumbers 
kd/n = 0.38, 0.57 and 0.76. The experimental data (23 
(see blue line (left axis) within the inset of Fig. [2] [331 ]) 
did indeed seem to confirm the theoretical predictions. 
However, various questions still remain to be clarified - 
the most prominent being the nature of the correspond- 
ing eigenstates and how they are structurally affected by 
these potential correlations. 
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FIG. 2: (Color online) Scaled inverse delay times for mi- 
crowaves propagating in ID waveguide. The different symbols 
correspond to energies \E{k = 0.57r)| < 1, \E(k = 0.7ir)| > 1 
being on both sides of the critical wavevector k = 0.57-7T. A 
nice data collapse is observed, indicating that in both cases, 
the statistical properties of delay times (and thus the struc- 
tural properties of wavefunctions) are unaffected by the cor- 
relation and correspond to exponentially localized wavefunc- 
tions, albeit the localization length for k = 0.5-7T is much larger 
than for k = 0.7n. This is reflected in the overall scaling pa- 
rameter {t^ 1 }. The dashed line is the result of the best fit 
from Eqs. (|12I13[1 , Inset: The experimental transmission co- 
efficient showing pass and stop bands is displayed by the blue 
line (left axis). The values for (Tj~ f ) are shown by the red 
circles (right axis). [33| . 



With the help of the iteration relations, Eqs. (|5l6|) . we 
have investigated the scaling properties of (t^ 1 ) for the 
correlated model, Eqs. (|7l9[ ). Two energies E from both 
sides of the mobility edge k = 0.577T have been cho- 
sen. In Fig. [2] we report our numerical data by re- 
ferring to the scaling variables /3_i and A_i, defined 
in Eq. ([T]). The data correspond to various system 
sizes L e 10\ 10 2 , 10 6 , 10 7 and disordered strengths 
e G {0.1, 0.5, 2.5, 5}. The remarkable agreement between 
the data from both sides of the "mobility" edge con- 
firms again the theoretical prediction, Eq.([T]), and indi- 
cates clearly that the corresponding cigenfunctions have 
the same structural properties, thus being unaffected by 
the potential correlations. Using the scaling properties 
of the Wigner delay times, we are able to conclude that 
k = 0.57-7T does not correspond to any mobility edge sepa- 
rating extended from exponentially localized eigenstates. 
Rather in both energy regimes the eigenstates are struc- 
turally the same (i.e. exponentially localized), albeit 
the localization length is drastically different. This is 
reflected in the overall scaling factor (t^ 1 ) (used to scale 
the data according to Eq. {T])), illustrated by the red 
circles (right axis) within the inset of Fig. [2l Note that 
Too ~ ^oo (see for example Q). As we can see from Fig. 
[21 at the pass-band region, (t^ 1 ) is much smaller than 
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that of the stop-band region; i.e, is much larger in 
the former case, but nonetheless remains finite (a "true" 
transition would imply that (t^ 1 ) ~ L™ 1 and thus by 
increasing the system size the scaling factor had to go 
to zero). This abrupt change in the magnitute of (t^ 1 ) 
arround k ~ 0.57-7T is a fingerprint of the correlations 
imposed to the disordered potential. Nevertheless, after 
rescaling the data the universal scaling law Eq. (fTJ) is 
again satisfied. 

Disordered Optical Lattices - It was recently proposed 
in 34[ that we can observe Anderson localization of 
ultra-cold atoms scattered off a gas of atoms of another 
species or internal state, randomly trapped at the nodes 
of an optical lattice. Within this set-up, cooled vibra- 
tional ground-state atoms trapped at the nodes of a pe- 
riodic optical lattice act as (static) delta scatterers pro- 
vided that the kinetic energy of the incoming particles 
is less than the vibrational energy of the trapped scat- 



terers i.e. 



2m m 



The mathematical 



model that describes the motion of the incoming particle 
along the lattice direction is the Kronig-Pcnney model, 
Eqs. (|7l8p . in this case with binary on-site potential distri- 
bution. Localization is then dependent on three param- 
eters: wavevector k, disorder strength e, and the filling 
factor p G [0, 1]. The latter dictates a binomial distribu- 
tion of the on-site potential 



e, Cn < P 
0, Cn > P 



(11) 



where £„ is a random number given by a uniform dis- 
tribution and e is the disorder strength (34|. In the 
numerical simulations presented in Fig. [31 we used dis- 
order strengths e G {4.556, 0.5} and filling factors p G 
{0.01,0.025,0.05,0.1,0.9}. The larger disorder strength 
corresponds to numerical values used in [3~j |. The very 
nice overlap of the scaled delay times are once more in 
excellent agreement with the universality of the scaling 

law, Eqs.CEJE]). 

Universal Behavior - It is illuminating to plot all our 
numerical data in the new variables 



Y- q = In 



, X_ g =ln(A_J 



(12) 



In these variables the scaling for q = 1 has an extremely 
simple form 



Y-! = o_i + 



(13) 



with a-i ~ and b-i w 1. The data for the scaling 
in variables Y—i,X—i are presented in Fig. 0J The re- 
markable result is that the above simple scaling relation 
holds in a very large region of the scaling parameter, 
AX_i w 14. In fact, Eg. ([TBI is exact only for q = 1, 
corresponding to q' = 2 |l5l. Hij|. However, for other 
values of g, Eo. (|13p is still a good approximation (see in- 
set of Fig. 2] for the case q = 2). Placing Eq. ([T2"|) into 
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FIG. 3: Scaled inverse delay time for the Disordered Optical 
Lattice System, with e G {4.556, 0.5} and filling factor p G 
{0.01, ,0.025,0.05, 0.1,0.9}. The nice data collapse confirms 
the universality of the scaling law, Eos. (|ll2 |l. The dashed line 
is the result of the best fit from Eqs. (|12ll3[l , 



Eq.{TnD, we find that 



A_i (sec dashed lines in 



1-/3-1 _ 

Figs. [Tl [2l3|) . This takes the form (|2|) once we substitute 
for /3_i and A-i the expressions in Eq.([T]). In the inset 
of the same figure we also report the q = 2 moment of 
delay times by making use of the variables of Eq. (fT2"|) . 
The nice data collapse reconfirms the validity of Eq. fTS"]) 
where again a_2 ~ and 6_2 ~ 1 (note however that for 
q = 2 small deviations from the straight line are evident 
around X^2 ~ 0) 
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FIG. 4: Scaling of Eq.([TT2|l in the variables from Eq.jrjJ. 
Inset: Same as in the main figure but now for the q = 2 case. 

In conclusion, we have investigated the scaling proper- 
ties of inverse moments of Wigner delay times. We have 
shown that they are dictated by the scaling law, Eq.([T]), 
which can be rewritten in a more familiar way, Eq.@, re- 
sembling the scaling relation for the information lengths 



5 



of wavcfunction components. Our theoretical arguments 
have been tested in various physical models where the ap- 
plicability of the non-linear a model is questionable, thus 
strongly supporting the relation between wavefunction 
moments and inverse moments of Wigner delay times. 
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